Optimizing thermoelectric properties of filled 
MPt 4 Ge 12 _ :r Sb :r skutterudites by band engineering 
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On the basis of density functional theory (DFT) calculations thermoelectric properties are derived 
for Sb-doped skutterudites MPt4Gei2- :E Sb a ; with M=Ba,La,Th. It is predicted that the originally 
very small absolute values of Seebeck coefficients \S\ of the undoped compounds is increased by 
factors of 10 or more for suitable dopings. The optimal dopings correspond to a "magic" valence 
electron number for which all electronic states up to a (pseudo)gap are filled. The theoretical findings 
are corroborated by measurements of S for LaPt4Gei2-a ! Sb a ; skutterudites. DFT derived vibrational 
rattling-like modes for LaPt4Gei2 indicate a small value for the lattice thermal conductivity which 
in combination with a large value of S 2 makes the La-based skutterudites appear as promising 
thermoelectric materials. 

PACS numbers: 71.20.-b, 72.15.Jf, 72.20.Pa 
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Thermoelectric materials are of scientific as well as 
technological interest in particular because of their po- 
tential for environment friendly electric power genera- 
tion. The thermoelectric performance of a material is 
characterized by the dimensionless figure of merit, 

ZT = S 2 Ta/(n e i + n ph ), (1) 

where S represents the Seebeck coefficient, a the elec- 
trical conductivity, and K e i , n p h are the thermal conduc- 
tivities attributed to electrons and phonons, respectively. 
Skutterudites appear as promising thermoelectric mate- 
rials by combining suitable electronic, vibrational and 
thermal properties. This is due to their crystal structure 
which contains large voids formed by framework atoms. 
Filling the voids by heavier elements may reduce n p h due 
to the appearance of low-lying optical (so-called rattling) 
modes [l|. Additionally these filler atoms may provide 
an appropriate number of valence electrons for making 
S 2 large. At present, a ZT at elevated temperatures of 
about 1.6 is achieved by filled pnictide skutterudites 
Very recently, values of ZT = 1.7 at 800 K were reported 
for the unfillled skutterudite Fe2.sCoi.2Sbi2 under severe 
plastic deformation Q. 

Ge-based skutterudites MPt 4 Gei2 were recently syn- 
thesized [3|-[8j for which the cage framework is formed 
entirely by Ge atoms. Their Seebeck coefficients are 
rather low as compared to pnictogen-based skutter- 
duites [!, H, ®- 

The aim of the present work is to utilize density func- 
tional theory results for understanding the physical prop- 
erties of Ge-based and Sb-doped filled skutterudites in 
order to provide a recipe for improving their thermoelec- 
tric performance. The quantities provided by DFT are 



the electronic structure, the Seebeck coefficient, the en- 
ergy of formation as a function of doping and vibrational 
properties when searching for rattling modes. 

Previously, efforts have been made for investigating the 
electronic structure in order to optimize the thermoelec- 
tric performance. Mahan and Sofo [lj| proposed that 
a delta-shaped transport distribution function originat- 
ing from very narrow electronic density of states (DOS) 
around the Fermi energy Ep is required to maximize ZT, 



as also recently hinted by experiment [111 ]. For this pur- 
pose, doping is essential, as it will be applied also in 
the present study. On the other hand, the required sin- 
gularity of the DOS is not accomplishable in many real 
thermoelectric materials although they possess promis- 
ing thermoelectric properties [12H15I ] . Based on a tight- 
binding concept it was argued by Zhou et al. [l6| that ex- 
tremely narrow bands are not optimal for maximizing ZT 
but rather the type of scattering mechanism of the charge 
carriers is of importance. Here, we propose a scheme of 
engineering the electronic structure of MPt4Gei2 com- 
pounds by investigating the substitutionally doped com- 
pounds MPt4Gei2- a; Sb 2 ;. In fact, we predict an increase 
of the Seebeck coefficient by a factor 10 and more upon 
doping. These findings are corroborated by new experi- 
mental results for LaPt4Gei2- a; Sb a ;. 

According to Boltzmann's transport theory, the See- 
beck tensor is defined as 



S(T, N) 



1 



(2) 



with the electronic charge — \e\. The total number of 
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electrons N defines the chemical potential n by 

/oo 
n(E)f(E,fi,T)dE = N (3) 

when integrating the DOS n(E) weighted by the Fermi- 
Dirac function / = 1 / (1 + exp((E — /i) / kT) . In the spirit 
of a rigid band model (i.e. n(E) remains unchanged upon 
varying N) the components of the tensor S become now 
functions of N and temperature T. The tensor K n is 
associated with the electronic band structure by 



K n = --^r V V; (k) <g) V, (k)n (k) X 

477 (4) 

( £i (k)-^r/'( £i (k)) 



in which £i(k) and Vi(k) are the energy eigenvalue and 
the corresponding band velocity for vector k and for band 
i. The quantity T^(k) represents the corresponding relax- 
ation time comprising all scattering events. The function 
f'(E) = df jdE is the derivative of the Fermi-Dirac dis- 
tribution. The Seebeck tensor is diagonal with all com- 
ponents equal if the system has cubic symmetry or is 
symmetrically averaged, which we assume from now on. 
In fact, skutterudites crystallise in a cubic crystal struc- 
ture. Within Mott's approximation derived for low tem- 
peratures 



S 



17j the Seebeck coefficient is expressed as 
dn(E) 
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v 2 {E) dE t(E) dE 
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relating S to the DOS by the first term on the right side. 
This term becomes large for a large slope and a small 
value of n(Ep) which could be the case e.g. for doped 
semiconductors. Assuming a constant relaxation time t 
(i.e. the third term is zero) and a parabolic behaviour of 
the valence bands according to e(k) = E v — ti 2 k 2 /2m and 
of the conduction bands according to e(k) = H 2 k 2 /2m — 
E c one arrives at the relations 

k B T 



S n = — I const I 



Ef — E c 



for an n-type semiconductor and 

k B T 



S u 



\const\ 



E v — Ef ' 



Ef > E a 



Ep < E v 



(6) 



(7) 



for an p-type semiconductor. There, E c and E v are the 
bottom of the conduction band and the top of the valence 
band, respectively. Actually, a parabolic band behavior is 
observed for MPt4Gei2 (M=Ba,La,Th) skutterudites be- 
low the distinctive (pseudo)gaps, as revealed by Figs. [2] 
Maximizing IS"! in the spirit of Mott's approximation re- 
quires that Ep is close to the conduction or valence band 
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FIG. 1: DFT derived Seebeck coefficients S for MPt 4 Gei 2 , 
(M=Th,La,Ba). Left panel: S as a function of electron doping 
AN at T = 300 K; right panel: S as a function of temperature 
for the undoped compounds (i.e. AiV = 0). Doping AiV is 
defined with respect to MPt4Gei2, and AiV = 4,5,6 places 
the Fermi energy of MPt4Gei2 in the (pseudo)gap. S for 
AiV 7^ calculated within the rigid band model. 



edge and that the electron velocities are large according 
to the second term on the right side of Eq. [SJ 

For calculating the Seebeck coefficient according to 
Eqs|5] and HI within the constant relaxation time approxi- 



mation a modified version of the BoltzTrap program 18 1 
is used. Density functional theory (DFT) calculations 
were carried out by ap plying t he Vienna Ab initio Simu- 
lation Package (VASP) [19ll20j . The exchange correlation 
functional is approximated by the generalized gradient 
approximation as parametrized by Perdew, Burke and 
Ernzerhof [2l[ , and the pseudopotentials are constructed 
by Blochl's projector augmented wave method [22 , 23 1 . In 
all calculations the structures were fully relaxed utilizing 
a 5 x 5 x 5 Monkhorst and Pack k-point grid [24| ensur- 
ing accurate results. After structural optimization the 
eigenvalues £j(k) were derived on a very fine 25 x 25 x 25 
k-mesh as needed by BoltzTrap. 

Substituting Ge by Sb leads to a large homogeinity 
region in the phase diagram, and because of that it is 
unfeasible to perform fully relaxed DFT calculations for 
all possible structures. Therefore, such calculations were 
only made for the most important structures (i.e. the 
ones with the highest number of symmetry equivalent 
atomic configurations) , and out of this selection the one 
with lowest total energy was chosen for finally deriving 
the thermoelectric properties. 

The DFT derived Seebeck coefficients for MPt 4 Gei 2 
(M=Ba,La,Th) in Fig. [T] are rather small with values less 
than 10 (J.V/K being consistent with experiment 0, 
The left panel of Fig. Q] showing S as a function of elec- 
tron doping reveals that \S\ increase strongly when the 
Fermi level falls into the gap of the electronic structure, 
corresponding to a "magic" doping of AN m = 4,5,6 
for MPt4Gei2 (M=Th,La,Ba) (see also the corresponding 
panels in Fig. [2]). In particular the sign of S and its be- 
havior around the gap (i.e. change of sign and near diver- 
gency) are reasonably described by the low-temperature 
approximation for S of a p- and n-conducting semicon- 
ductor according to Eqs. [S] and [7J In Fig. [5Ja) the 
electronic structure at and below the (pseudo)gap re- 
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FIG. 2: Electronic structure of Pt4Gei2 (top left), BaPt4Gei2 
(top right), LaPt4Gei2 (bottom left), and LaPt4Ge7Sbs (bot- 
tom right). 



SmPt4Gei2 by previous calculations [25| . 

Another way of doping and adding valence electrons 
could be chosen by substituting the transition element. 
Noble metal elements such as Cu,Ag,Au are possible can- 
didates providing one more electron than the latest tran- 
sition elements Ni, Pd, Pt. The question is, how the 
electronic structure is affected and if the gap feature is 
preserved. When testing the replacement of Pt by Au 
by studying MAu4Gei2 compounds it turns out that the 
lowering of the 5d-level of Au in comparison to the Pt 
5d-levels destroys the (pseudo)gap feature and hence also 
the promising thermoelectrical properties. 
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veals a band with strong dispersion with its maximum 
in T (around 1.5 eV for the MPt4Gei2 compounds, and 
at Ef for LaPt4Sb5Gey), which behaves rather linearly 
along the T - H direction. These features resemble those 
of CoSb3, which is a narrow-gap semiconductor. As- 
suming that substitution of Co and Sb by similar ele- 
ments such as Pt and Ge does not significantly change 
the (pseudo)gap formation properties a simple electron 
counting rule can be constructed. Co 4 Sbi2 has 96 valence 
electrons whereas Pt4Gei2, BaPt4Gei2, LaPt4Gei2 and 
ThPt4Gei2 have only 88, 90,91 and 92 valence electrons, 
respectively. For placing now Ep into the (pseudo)gap 
the Pt-Ge compounds have to be doped by 8,6,5 and 4 
electrons ("magic" doping), correspondingly. 

Figure [2] reveals rather minor influences on the 
(pseudo)gap when filling the voids of Pt4Gei2: it is pre- 
served for Ba- and LaPt4Gei2 and becomes a true but 
small gap for ThPt4Gei2 (not shown). For LaPt4Gei2 
and ThPt4Gei2 bands above the conduction band edge 
are rather flat due to their localized f-like character. 
Then, St q^ is becoming large and as indicated by Eq. [5] 
IS | increases accordingly. Apart from changing the sign, 
this is the reason for the asymmetry of S below and above 
the (pseudo)gap, with 15*1 significantly larger for dopings 
A A > AN m . The crucial point is to find appropriate 
filler atoms which provide the correct number of doping 
electrons for Ep falling into the (pseudo)gap. 

When the filler atom M is chosen for doping, an ele- 
ment which provides 8 valence electrons is needed. Then, 
the dopant must be a lanthanide or actinide atom with 
f-states as valence states. However, due to their local- 
ized nature the f-states may strongly influence the elec- 
tronic structure at Fermi energy. For example, a metall- 
lic rather than semiconducting behavior is predicted for 



FIG. 3: DFT derived Seebeck coefficients as a function of 
doping at T=300K (left panel) and temperature (right panel) 
for selected MPt4Sbi2- a: Ge a; compounds. The change of sign 
of S occurs at the "magical" composition, for which the total 
number of valence electrons is 96. 
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FIG. 4: DFT derived (left panel) and experimental Seebeck 
coefficients (right panel) as a function of temperature for 
LaPt4Gei2- a; Sba ; compounds. 

A third concept for substitutions is the replacement 
of Ge by elements such as Sb with a larger number of 
valence electrons than 4. Filled skutterudites such as 
BaPt 4 Sb 6 Ge 6 , LaPt 4 Sb 5 Ge 7 and ThPt 4 Sb 4 Ge 8 have the 
"magic" number of valence electrons, namely 96, and 
what is important for them the desired electronic struc- 
ture properties are preserved, as shown in Fig. [2] for the 
La compound as an example. Even a small gap opens 
up, which is not the case for undoped LaPt 4 Gei2- For 
the calculated thermoelectric properties of LaPt 4 Sb 5 Ge 7 
it turns out that \S\ is substantially increased by more 
than a factor of 10 (Fig. [3]) reaching now values even 
larger than 100 /iV/K for the La-compound and more 
than 200 /iV/K for the Th-compound. Figure Mb) de- 
picts temperature dependent Seebeck coefficients of some 
MPt4Gei2-a;Sb x compounds around the critical compo- 



4 



sition (i.e. "magic" doping). Values of |5| larger than 100 
/iV/K are reached in a temperature range of 400-800 K 
depending on the compound and its composition. From 
these findings one may speculate that in general large 
Sccbcck coefficients are obtained with a proper doping of 
Sb for other Ge-based skutterudites. 

Figure [4] compares DFT derived and experimental re- 
sults for LaPt4Gei2-xSb2; compounds. Clearly, the ex- 
perimental evidence corroborates the theoretical predic- 
tion of a very strong increase of |5| upon doping by Sb. 
This enhancement effect seems to be more pronounced 
by the DFT results, which predict a maximum of S(T) 
ranging from 500 to 800 K for p-type materials depend- 
ing on the amount of substitutions. The DFT data also 
show the sensitivity upon substitution. Increasing the 
amount of Sb above the critical value (i.e. crossing the 
gap) changes its sign and results in comparable (although 
slightly smaller) values of \S\. 

For maximizing the figure of merit in Eq. [1] the power 
factor S 2 a should be large. Indeed, from the DFT calcu- 
lations it was found that for LaPt4Ge7.iSb6.g the power 
factor reaches a maximum at 600K, at a temperature 
which is useful also for technological applications. It 
should be noted that an absolute value of S 2 a cannot 
be obtained from the first-principles calculations, be- 
cause the relaxation time (assumed to be constant) is 
not known. For the Seebeck coefficient, however, the 
constant relaxation time cancels out. Nevertheless, an 
enhancement of at least two orders of magnitude in the 
power factor is achieved for each compound within the 
constant relaxation time approximation. Experimental 
data for LaPt4GeySb5 at 800 K reveal a power factor of 
S 2 p = 7.3* 10- 5 W/K 2 . 

According to Eq. [T] a large figure of merit can also 
be achieved when the thermal conductivities are small. 
Concerning n e i, the electronic contribution, it should be 
proportional to the DOS at Fermi energy. Inspecting 
Fig.[2]for LaPt 4 Ge 7 Sb 5 it is obvious that the DOS is very 
small just below the gap. Therefore, experimenting with 
the Sb substitutions close to the critical Sbs composition 
might lead to a very small value of n e i. Concomitantly 
the electrical conductivity a becomes also small, which 
counteracts the enhancement effect of n e i for the figure 
of merit. 

Considering the phonon contribution to the thermal 
conductivity the La-compound seems very promising to 
yield a low value of K P h, because for LaPt4Gei2 well local- 
ized rattling-like vibrational modes were found by DFT 
calculations which resulted in a pronounced peak of La- 
character at 1.7 THz in the phonon DOS. Replacing La 
by Ba, which has a very similar mass but one valence 
electron less, results in a very different phonon spectrum 
at lower frequencies without any hint of a rattling-type 
feature. Results of these calculations of vibrational prop- 
erties of MPt4Gei2 compounds will be presented in a 
forthcoming work [26J. 



For synthesizing the material it is important to esti- 
mate its thermodynamic stability. Doing this by deriving 
the DFT enthalpy of formation AH it is found that about 
1 kJ/mole is gained per substituting Sb atom with the 
minimum of AH (i.e. strongest bonding) for the doped 
compound MPt4Ge7Sbs. Further increasing the amount 
of Sb reduces the bonding and hence the energy of for- 
mation is becoming less negative. This behavior reflects 
the filling of bonding valence states upon doping until 
the gap is reached for 5 Sb atoms (Fig. [21 right bottom 
panel). Adding more Sb adds non- or antibonding con- 
duction band states to the total energy of the compound 
and hence AH is getting less negative. 

On the basis of DFT calculations we predict new and 
promising thermoelectric materials, when in MPt4Gei2 
skutterudites Ge is suitably replaced by Sb. In particu- 
lar, LaPt4Ge7_ x Sb5_ x is studied which reveals strongly 
enhanced values of the Seebeck coefficient up to about 
200 fx V/ K at 600 K. The "magic" composition of 7 
Ge and 5 Sb atoms can be understood in terms of band 
engineering revealing the importance of electronic bands 
with large dispersions below a (pseudo)gap, very simi- 
lar to CoSb3. Samples of LaPt4Gc7_ :E Sb5_ x were ex- 
perimentally synthesized. Measurements of the Seebeck 
coefficient strongly supports the DFT predictions. In 
addition, from the DFT derived vibrational properties 
of LaPt4GeySb5, which reveal rattling modes, it can be 
concluded that the phonon thermal conductivity is ex- 
pected to be small, which -in addition to the large See- 
beck coefficient- would enhance the figure of merit. 
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